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ABSTRACT 



Context, The rest-frame near-IR spectra of intermediate age (1-2 Gyr) stellar populations are dominated by carbon based absorption 
features offering a wealth of information. Yet, spectral libraries that include the near-IR wavelength range do not sample a sufficiently 
broad range of ages and metallicities to allow for accurate calibration of stellar population models and thus the interpretation of the 
observations. 

Aims. In this paper we investigate the integrated J- and //-band spectra of six intermediate age and old globular clusters in the Large 
Magellanic Cloud (LMC). 

Methods. The observations for six clusters were obtained with the SINFONI integral field spectrograph at the ESO VLT Yepun 
telescope, covering the J (1.09 - 1.41 yum) and //-band (1.43 - 1.86 /mi) spectral range. The spectral resolution is 6.7 A in J and 6.6 A 
in //-band (FWHM). The observations were made in natural seeing, covering the central 24" x 24" of each cluster and in addition 
sampling the brightest eight red giant branch (RGB) and asymptotic giant branch (AGB) star candidates within the clusters' tidal 
radii. Targeted clusters cover the ages of -1.3 Gyr (NGC 1806, NGC 2162), 2 Gyr (NGC 2173) and ~ 13 Gyr (NGC 1754, NGC 2005, 
NGC2019). 

Results. //-band C2 and A"-band 12 CO (2-0) feature strengths for the LMC globular clusters are compared to the models of Maraston 
(2005). C2 is reasonably well reproduced by the models at all ages, while 12 CO (2-0) shows good agreement for older (age > 2 Gyr) 
populations, but the younger (1.3 Gyr) globular clusters do not follow the models. We argue that this is due to the fact that the empirical 
calibration of the models relies on only a few Milky Way carbon star spectra, which show different 12 CO (2-0) index strengths than 
the LMC stars. The C2 absorption feature strength correlates strongly with age. It is present essentially only in populations that have 
1-2 Gyr old stars, while its value is consistent with zero for older populations. The distinct spectral energy distribution observed for 
the intermediate age globular clusters in the /- and //-bands agrees well with the model predictions of Maraston for the contribution 
from the thermally pulsing asymptotic giant branch phase. 

Conclusions. In this pilot project we present an empirical library of six LMC globular cluster integrated near-IR spectra that are 
useful for testing stellar population models in this wavelength regime. We show that the //-band C2 absorption feature and the J-, 
//-band spectral shape can be used as an age indicator for intermediate age stellar populations in integrated spectra of star clusters 
and galaxies. 
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1. Introduction 

With the appearance of more and better near-IR instruments, the 
studies of the stellar populations in galaxies in different environ- 
ments and at various redshifts, exploring the rest-frame wave- 
leng th region between 1 and 2.5 yum, have become more frequent 
(e.g. ISilva et al.ll2008l lEancon et al.ll200"8t iHvvonen et al.11200 1 



Marmol- Oueralto et all 120091: ICesetti et all 120091: Martins et al 
2010t iRiffel et al.ll201 lah~ Thus a reliable framework to inter- 



pret these observations is needed. The comparison of observa- 
tions with stellar population models is an approach that proved 
to be successful in the optical wavelength regime over the last 
decades. However, the near-IR range of the evolutionary pop- 
ulations synthesis (EPS) models still lacks a proper empirical 
calibration. Recent advances in theoretical as well as empiri- 



* Based on observation collected at the ESO Paranal La Silla 
Observatory, Chile, Prog. ID 078.B-0205 



cal calibration include in particular work on asymptotic giant 
branch (AGB) evolutionary phase, whose contribution is cru- 
cial to the near-IR part of the models, (e.g . iMarigo etaL 20081: 
iGirardi et al.l 120 1 Ok iMarastonl 1 1 998L 120051) . In iLvubenova et all 



(1201 (A hereafter Paper I) we described our efforts to create a 
near-IR library of integrated spectra of globular clusters to serve 
as a cali brator and test bench for existing and future EPS mod- 
els (e . g. iBruzual & CharTo1l2003t lMarastonll2005l IConrov et all 
120091: iMaraston & Strombackll2011l) . For this pilot project we 
chose six globular clusters in the LMC, for which detailed infor- 
mation about the age and chemical composition exists, based on 
their resolved light. Our sample consists of three old (> 10 Gyr) 
and metal poor ([Fe/H] —1.4) clusters, namely NGC 1754, 
NGC 2005, and NGC 2019, as well as three intermediate age 
(1 - 2 Gyr) and more metal rich ([Fe/H] ^ -0.4) clusters: 
NGC 1806, NGC 2162, and NGC 2173 (see TableQ}. 
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Table 1. Target globular clusters in the LMC. 



Name 


Ape (Crvri 


IFc/HI 


'-mosaic v^aJ/ 


1 II 


1 1 


S/N; 


S/Nw 

.J/1 >tl 


(i) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


NGC 1806 


1.1' , 1.7 s , 1.5 s , 1.9 i ,0.5 p 


-0.23*, -0.7 V 


3.4 x 10 4 






>60* 


>45* 


NGC2162 


l.l 1 , 1.3 3 ,2 2 ,0.9\ 1.25 5 


-0.23* , -O^ 


0.5 x 10 4 


21'.' 37 


19772 


>45* 


>40* 


NGC 2173 


2 1 , 2.1 3 , 1.6 6 , 4.1 2 , 1.5 7 , l.l 4 , 1.6 5 


-0.24* , -0A2 f , -0.51* 


0.9 x 10 4 


34'.' 35 


39375 


>80* 


>45* 


NGC 2019 


10 1 , 16'°, 13.3 11 , 16.3°, 17.8* 


-1.23", -1.18*, -1.37', -1.10'' 


6.7 x 10 4 


9772 


12176 


95 


150 


NGC 2005 


10',6.3 10 , 16", 15.5", 16.6* 


-1.35", -1.92*, -1.80 c , -1.33 d 


3.9 x 10 4 


8765 


9878 


140 




NGC 1754 


10 1 , 7'°, 14" , 15.6"-* 


-1.42", -1.54* 


2.7 x 10 4 


1172 


14279 


80 





N otes. (1) Cluster name, (2) age of the cluster in Gyr, derived usin g diffe rent methods: (1> iFrogel et all d 19901) - based o n th e SWB type; 
(2) iLeonardi & Rosd d2003h - i ntegrated spectroscop y; < 3) iGeisler et al.l lll997h . (4) iGirardi et al.l (Il995h. (5) [Kerber et all j2007h . (6) iBertelli et all 
( 120031) , (7) l Woo et al.lHo03h. m |Mackev et alj d2008h . <9) lMilone et all d2009l) - all CMDs. ""' iBeaslev et aljj2002l) - integrated spectro scopy (H/3- 
Mgb), (11) .Beas lev et alj d2002l) - in tegrated spectroscopy (Hy-(Fe)), (3) [Fe/H] der i ved u sing different methods: M lOlsen et al J d 1998f) - slope of 
the RGB, (b) lOlszewski et al J d 1 99lh - l ow-resolution Call triplet, (c) Ijohnson et all d2006|) - h igh-resolution Fel, {d) I Johnson et all d2006 ] ) - high- 
resolu tion Fe II. (e) iDrrschet al.1 d2000j) - Stromgren photometry, (fi lGrfx!r!olski"etal.l d2006l) - low-resolution Call triplet. w llVIucciarelli et all 
(2008) - high-resolution spectroscopy (4) sampled bolometric luminosity within the clusters centra l mosaics in Lq, (5) half-light radius and (6) 
tidal radius of the King-model cluster fit from the catalogue of McLaughlin & van der Marel (2005), (7) and (8) S/N of the integrated spectra in 
J- and //-band, respectively; (,) lower threshold, see Sect. l2.3l for an explanation. 



In Paper I, using the integrated /<f-band spectra of the globu- 
lar clusters, we demonstrated the feasibility of our observational 
approach, as well as discussed some discrepancies that arise be- 
tween recent EPS models and the data. We argued that the main 
reason for these discrepancies is the incompleteness of the cur- 
rent stellar spectral libraries of asymptotic giant branch (AGB) 
stars that are used to create the models. Also, we illustrated how 
the presence of luminous carbon-rich AGB stars significantly 
changes the observed spectrophotometric properties in the K- 
band, especially the 12 CO (2-0) 2.29 //m feature as predicted by 
iMarastonl (120051) . despite the observational discrepancies. 

This paper complements Paper I by presenting results based 
on the integrated J- and //-band spectra of the same globu- 
lar clusters. To date there are very few studies of the spectral 
properties of glob ular clusters in this wavelength regime (e.g. 
iRiffel etal.ll2011bh . The paper is organised as follows: in Sect. [2] 
we briefly recall our observing strategy and data reduction meth- 
ods. In Sect. [3] we compare the observed overall spectral energy 
distributions of the globular clusters to stellar population models. 
Sect. |4]is devoted to the C2 absorption feature in the //-band and 
its behaviour in the spectra of globular clusters. In Sect. |5]we dis- 
cuss potential reasons for the disagreement between model pre- 
dictions and our observations, based on stellar libraries. Finally, 
in Sect. [6] we present our concluding remarks. 

2. Observations and data reduction 

2.1. Observations 

Our sample of globular clusters, whose basic pr operties are listed 
in TableCD were observed w ith VLT/SINFONI ( lEisenhauer et al.l 
12003b | Bonn et et al ] I2004 in poor seeing conditions without 
adaptive optics in the period October - December 2006 (Prog. 
ID 078.B-0205, PI: Kuntschner). The detailed descriptions of 
the observing strategy, target selection and data reduction are 
given in Paper I, where we presented the /iT-band spectra from 
our project. In the present paper we explore the data obtained in 
the /- and //-bands, where minor differences in the data reduc- 
tion exist. Briefly, in order to sample at least the half-light radius 
of the clusters, we observed a 3 x 3 mosaic of the largest field- 
of-view that SINFONI offers (8" x 8"). Thus, the resulting cov- 
erage was 24" x 24", with a spatial sampling of 07 125 per pixel. 
For the observations we used the standard near-IR nodding tech- 
nique. We observed each mosaic tile three times with individual 



integration times of 50 s and dithering of 0725 to ensure better 
bad pixel rejection. We took sky exposures between each mosaic 
tile. Finally, within each cluster's tidal radius, we selected up to 
eight bright stars with colours and magnitudes typical for the red 
giant branch (RGB) and AGB stars in the LMC (see Table 3 in 
Paper I), which were observed in addition to the cluster mosaics 
in order to better sample the short lived AGB phase. The integra- 
tion times for these stars were 10 s. After each science target and 
at a similar airmass, we observed an A-B dwarf star for telluric 
absorption correction (Table 4 in Paper I). 

2.2. Data reduction 

We used the SINFONI pipeline version 2.0.5 to perform the ba- 
sic data reduction on the three exposures per mosaic tile plus 
two bracketing sky exposures. In brief, the pipeline extracts the 
raw data, applies distortion, bad pixel and flat-field corrections, 
wavelength calibration, and stores the combined sky-subtracted 
spectra in a 3 -dimensional da ta cube. Then on ea ch resulting data 
cube we ran the lac 3d code dDavies et alj2010l) . whose purpose 
is to detect and correct residual bad pixels that are identified us- 
ing a 3D Laplacian edge detection method. This code not only 
removes residual bad pixels, but also produces an error data cube 
for each input science data cube, which is very useful, because 
the SINFONI pipeline does not provide an estimate of the error 
propagation during data reduction. The derived signal-to-noise 
ratio, using these error spectra, is in agreement with the signal- 
to-noise ratio deriv ed using an empirical method described by 
IStoehret al.l(l2007h . 

In the /-band we had to switch off the pipeline option 
to correct for sky line r esiduals based on the algorithm de- 
scribed in iDaviesi (120071 i.e. --objnod-sky_cor = FALSE). 
This was needed due to the complicated pattern of overlap- 
ping OH and O2 sky emission lines in this wavelength region 
(see e.g. iRousselot et al.l l2~000h . which leads to their over- or 
under-subtraction when applying the sky residual correction. 
Therefore, we applied simple sky subtraction and indicated the 
regions where the sky lines are the strongest in Fig. [T](i.e. 1.21 — 
1 .236 //m and 1 .26 - 1 .295 //m). 

We reduced the telluric star data in the same way as the 
science frames. Then for each telluric star we extracted a one- 
dimensional spectrum, removed the hydrogen absorption lines 
fitted with a Lorentzian profile, and divided the star spectrum 
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Fig. 1. Left panel: /-band, integrated spectra of the six LMC globular clusters. The shaded areas indicate the regions with contami- 
nation by strong OH and O2 sky line residuals. Right panel: //-band, integrated spectra of four of the LMC clusters. Due to unstable 
atmospheric conditions, the //-band spectra of NGC 1754 and NGC2005 are heavily contaminated by sky line residuals and thus 
exibit a very low S/N ratio a nd are not shown. Each spectrum is normalised to i t s med ian value. Line identifica tions are based on 
the stellar spectral atlases bv lLancon & Rocca-Volmerangek 1992). W allace et al.l (l2000l) and lRavner et al.l (120091) . 



by a black body spectrum with the same temperature as the star. 
The last step in preparing the telluric correction spectrum was 
to apply small shifts (<0.05 pixels) and scalings to minimise the 
residuals of the telluric features. Finally, we divided each science 
data cube by the corresponding telluric spectrum. In this way we 
also obtained an approximate relative flux calibration. An ab- 
solute flux calibration was not possible due to non-photometric 
conditions. 

The telluric star HD 44533, used for the telluric correction of 
the observations for NGC 2019 and its surrounding stars, had an 
unusual shape of the hydrogen lines. It seems that this star shows 
some hydrogen emission together with absorption. Thus, to re- 
move the hydrogen lines, we interpolated linearly over the af- 
fected regions, which are not overlapping with any indices used 
in this study. However, we add a word of caution to treat these 
regions with care in case the spectra of NGC 2019 are used in 
the future for other purposes. 

2.3. Spectra integration 

To obtain the final, integrated J- and //-band spectra for the six 
globular clusters, we estimated the noise level in each recon- 
structed image from the central mosaic data cubes. We assumed 
that this noise is due to residuals after the sky background cor- 
rection. Thus, we computed the median residual sky noise level 
and its standard deviation, after clipping all data points with in- 
tensities of more than 3<x. We then selected all spaxels, which 
have an intensity more than three times the standard deviation 



above the median residual sky noise level. We summed them and 
normalised the resulting spectrum to 1 s of exposure time. The 
same approach was used to obtain the spectra for the additional 
bright RGB and AGB stars. 

The final step in preparing the luminosity weighted, inte- 
grated spectra for the six globular clusters was to add the ad- 
ditionally observed bright stars outside the central mosaics (for 
a detailed discussion on cluster membership and which stars 
where included, see Sect. 4 in Paper I and Table |2). Fig. Q] 
shows the final, integrated J- and //-band spectra for the six 
LMC globular clusters. The spectra of the three intermediate 
age clusters are shown on the top (NGC 1806, NGC 2162, and 
NGC 2173). These spectra display numerous absorption features 
due to carbon based molecules, like CO, C2, and CN which 
makes them appear to be very noisy, while in fact they are 
not. The spectra of the old and metal poor clusters (NGC 2019, 
NGC 2005, and NGC 1754) are shown at the bottom. These 
are dominated by late-type giant stars, where the metallic lines 
are more prominent. To identify the spectral absorption lines 
displayed in Fig. [D we used the stellar spectral atlase s of 
[ Lancon & Rocca-Volmerangd(ll992h . lWallace et ail d2000l) , and 
iRavner et al. ( 20091) . The weather conditions during the observa- 
tions of the //-band spectra for NGC 1754 and NGC 2005 were 
not favourable, which resulted in a signal-to-noise ratio that was 
too low and contamination by multiple sky line residuals. Thus 
we decided to exclude these two spectra from our analysis. 

The spectral resolution of our spectra, as measured from 
arc lamp frames, is 6.7 A and 6.6 A (FWHM) in the / and 
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//-band, respectively. Finally, we estimated the S/N ratio of 
each integra t ed glo bular cluster spectrum using the method of 
IStoehr et al.l (120071) and listed the values in TableQ] This method 
allows us to compute the S/N from the spectrum itself. However, 
due to the numerous absorption features mimicking noise in the 
spectra of globular clusters containing carbon-rich AGB stars, 
i.e. the three intermediate age clusters marked with an asterisks 
in TableQ] we are able to give only a lower threshold. Based on 
the integration times and luminosities of the clusters, we con- 
clude that the quality of these intermediate age spectra is as 
good, if not better, than those of the old globular clusters. 

3. The overall near-IR spectral energy distribution 

Having completed the primary goal of our project, namely to 
provide an empirical library of integrated near-IR spectra of 
globular clusters, we show in Fig. [2] the overall J, H and K- 
band spectra as derived from our SINFONI observati ons, com- 
pared to the spectral energy distributions from the iMarastonl 
(120051) stellar population models. The observed spectral seg- 
ments were scaled individually relative to the model providing 
the closest match in age and metallicity and were normalised 
to unity at 1.25 pm (see Fig. caption). In the figure the mod- 
els are represented with solid, black lines. With the exception 
of the /-band part in old globular clusters, where we were un- 
able to achieve good sky background removal, as explained in 
Sec. 12.21 the models agree well with the general spectral shape 
of our spectra. Remarkably, the distinct spectral energy distribu- 
tion (the "sawtooth" pattern) in the /- and //-bands caused by 
the contribution from carbon-rich thermally pulsing asymptotic 
giant branch (TP-AGB) phase is visible in the spectral regions 
covered by SINFONI for the intermediate age clusters. As ex- 
pected the pattern is less prominent for the cluster NGC2173 
due to its slightly older age. Also, the most prominent features 
such as C2 at 1.77 pm and the CO bandhead at 2.29 pm are well 
matched. Thus, when they become available, it will be very in- 
teresting to make a comparison with higher spectral resolution 
EPS models. Weaker features, as the ones indicated in Fig. [Tj 
could be studied in detail. In this figure, in addition to carbon 
based molecular absorption bands, several other absorption fea- 
tures are indicated. They are mainly due to metallic lines, such 
as Fe I, Mg I, Si I, and Al I. However, due to their relative weak- 
ness and associated difficulty of measuring them in galaxies, we 
decided not to discuss them further in this paper. 
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Fig. 3. Top panel: Spectrum of NGC 1806 with overplotted cen- 
tral passband (red) and the clas sical continuum region (blue) for 
the C2 index ( IMarastonl 120051) . In orange we overplot the new 
continuum region for the C2 index defined in this paper (see 
Sect. H]i. Bottom panel: Comparison between the classical and 
the redefined C2 index, measured on our own LMC star obser- 
vations (solid black and red symbo ls), as well as the stars from 
the Milky Way spectral libraries of lLancon & Mouhcinel (120021 
orange ) and lRavner etail (120091 black diamonds, green squares, 
red asterisks). The dashed line shows the one-to-one relation, the 
solid line represents a linear least squares fit to all data points. 



4. The //-band C 2 index 

One of the most prominent and easy to quant ify features in the 
//-ba nd is the C 2 (0 - 0) bandhead at 1 .77 pm (iBallik & Ramsavl 
119631) . The existence of the Co molecule is typical for carbon- 
rich stars, where the ratio of carbon to oxygen (C/O) atoms 
is larger than 1. This type of star is the main contributor to 
the near-I R light of intermediate age (1-3 Gyr) stellar popula- 
tions (e.g.lFerraro et al.ll995tlGirardi & Bertellill998l IMarastonl 



I1998L 12005 ). Thus, it is expected that the C2 absorption feature 
will be strong in globular c lusters and galaxies with stellar popu - 
lations in this age interval ( lLancon et al.lll999tlMarastonll2005l) . 
Also, more metal poor stellar populations exhibit a stronger 
C2 index on average, be cause dredge-up more easily leads 
to C/Q>1 in that regime i Wagenhuber & Grpenewegenl 1 19981: 
iMouhcine & Lanconlfeoollweiss & Fergusonll2009l 

The C2 index reflects the de pth of the C2 absorption feature 
at 1.77jt/m dAlvarez et al. II2000I) . In our work we used the def- 
inition employed by IMarastonl (120051 hereafter we refer to this 



definition as "classical"). The index is defined in magnitudes, as 
the flux ratio between a central passband (1 .768 - 1 .782 pm) and 
a continuum passband (1 .752 - 1 .762 pm), and is normalised to 
Vega by subtracting a zero point of 0.038"'. In this definition of 
the index, the continuum band is on top of H2O and 12 CO fea- 
tures (see Fig. [3] top panel). While carbon-rich stars are not ex- 
pected to have H2O absorption, these two f eatures may vary w ith 
the pulsation period for oxygen-rich stars (tLoidl et al.ll200ll) . 

Because of these additional contributions to the continuum 
passband we explored a modified index definition by shifting 
the continuum passband to shorter wavelengths relative to the 
classical definition (1.74 - 1.75 pm, shown in orange in Fig. [3}. 
The Vega zero point to be subtracted is 0.037"'. We measured 
the two indices for the LMC stars from o ur own observations, 
the car bon-rich averaged stellar spectr a of Lancon & Mouhcinel 
(120021) and the stars from the library of lRavneret al.ld2009b . The 
comparison is shown in Fig. [3] (bottom panel). We found that for 
these samples the scatter is not significant, but one might expect 
it in a larger sample, and in particular if one includes AGB stars 
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Fig. 2. J, H, /f-band spectral energy distributions of our sample of six LM C globular cluste rs. The black solid lines indicate the 
most closely matching (in age and metallicity) stellar population model of Marastonl (I2005I) : i.e. age = 13Gyr, [Z/H] = -1.35 
for NGC 1754, NGC2005, and NGC2019, age = 2Gyr, [Z/H] = -0.33 for NGC2173, and age = 1.5 Gyr, [Z/H] = -0.33 for 
NGC 1806 and NGC2162. The spectra are normalised to the flux at 1.25 fim. Note the prominent spectral shape differences in the 
J- and //-bands between old and intermediate age clusters. The shaded areas indicate the regions with contamination by strong OH 
and O2 sky line residuals. 



Table 3. Measurements of Deo and C2 indices for the LMC 
globular clusters. 



Name 



D c 



C 2 (mag) 



NGC 1806 
NGC 2162 
NGC 2173 
NGC 2019 
NGC 2005 
NGC 1754 



1.129±0.005 
1.108±0.010 
1.186±0.005 
1.068±0.003 
1.086±0.003 
1.082±0.005 



0.214±0.010 
0.235±0.022 
0.068±0.012 
0.006+0.004 



Notes. The Dm index is defined in 



ilMarmol-Oueralto et all d2008l) and 
the C2 index is taken from lMarastor (l2005h . 



in different pulsation phases. For the present study we decided to 
use the classical index definition for direct comparison to mod- 
els. The measurements for the stars in our library are listed in 
Tableland for the globular clusters in Table|3] 



Fig.|4]shows a comparison between our measurements of the 
C2 index f or the LMC globu lar clusters and the SSP model pre- 
dictions of lMarastonl (l2005h . In the top panel, for guidance, we 
repeat Fig. 9 from Paper I, showing the /if -band Deo index. This 
index is sensitive to the metallicity at older ages, but also shows 
a dependence on the presence of carbon-rich stars at interme- 
diate ages. We measured the Deo an d C2 indices at the nomi- 
nal spectral resolution of the data. Current models have lower 
spectral resolution than our data. However, due to the broadness 
and intrinsic strength of these spectral features, the differences 
in spectral resolution will not affect our conclusions. We tested 
this by broadening our globular cluster spectra to spectral res- 
olutions ranging from 50 to 400 km s in steps of 50 km s . 
The largest offset we measured was -0.03 for both the Deo an d 
C2 indices, which is much smaller than the observed differences 
between observations and models. 

For old stellar populations (>3 Gyr) the models predict an 
approximately constant and zero value of the C2 index (Fig. [4] 
bottom panel). The C2 index measurement for the old globu- 
lar cluster NGC 2019 in this age range whose metallicity is best 
compared to the model with [Z/H] = -1.35 (solid, blue line), 
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Table 2. Near-IR colours and spectral indices of the additional bright RGB and AGB stars. 



Name 


K (mag) 


(J-K) 




Deo 


C 2 (mag) 


Notes 


Cluster 


(1) 


(2) 


(3) 




(4) 


(5) 


(6) 


(7) 


2MASS J04540 1 27-702634 1 


11.40 


1.13 


1 


.208+0.013 


0.034+0.010 




NGC 1754 


2MASS J0454077 1 -7024398 


11.86 


1.26 


1 


.223+0.021 


0.152+0.013 






2MASS J04543522-7027503 


12.26 


1.07 


1 


.143+0.015 


-0.011+0.018 






2MASS J04542864-7026142 


12.62 


0.91 


1 


.112+0.018 


0.052+0.025 






2MASS J04540935-7024566 


12.73 


0.96 


1 


.088+0.021 


-0.048+0.020 






2MASS J04541 1 88-7028201 


13.04 


0.91 


1 


.177+0.030 


0.045+0.009 




J) 


2MASS J04540536-7025202 


13.36 


0.93 


1 


,122±0.029 


0.156+0.012 






2MASS J0530222 1 -6946 1 24 


10.65 


1.21 


1 


.258+0.008 


-0.021+0.010 




NGC 2005 


2MASS J05300708-6945327 


11.66 


1.04 


1 


,213±0.012 


-0.044+0.018 






2MASS J05295466-6946014 


11.68 


1.12 


1 


.246+0.010 


-0.045+0.020 






2MASS J05300704-6944031 


11.80 


1.03 


1 


.251+0.010 


-0.039+0.027 






2MASS J05300360-694431 1 


12.02 


0.97 


1 


.216+0.013 


0.063+0.036 






2MASS J05295822-6944445 


12.05 


1.22 


1 


.233+0.016 


-0.040+0.041 






2MASS J05320670-7010248 


10.49 


1.13 


1 


.136+0.014 


-0.018+0.005 




NGC 2019 


2MASS J05313862-7010093 


10.56 


1.39 


1 


.196+0.017 


0.339+0.006 


C 




2MASS J05315232-7010083 


10.61 


0.94 


1 


.200+0.016 


-0.022+0.005 






2MASS J05321 152-7010535 


10.87 


1.06 


1 


.176+0.018 


-0.036+0.005 






2MASS J05321647-7008272 


10.95 


1.13 


1 


.183+0.016 


-0.008+0.006 






2MASS J0532041 8-7008 1 5 1 


11.22 


1.02 


1 


.204+0.024 


-0.047+0.007 






2MASS J05021232-6759369 


10.23 


1.84 


1 


.051 ±0.0 13 


0.458+0.007 


+,c 


NGC 1806 


2MASS J05020536-6800266 


10.63 


1.61 


1 


.139+0.015 


0.500+0.009 


+,c 




2MASS J05015896-6759387 


10.69 


1.76 


1 


.065+0.013 


0.481+0.010 


+,c 




2MASS J05021623-6759332 


11.02 


1.06 


1 


,211±0.012 


-0.047+0.007 


+ 




2MASS J05021870-6758552 


11.32 


1.11 


1 


.267+0.015 


-0.071+0.011 


+ 




2MASS J0502 1846-6759048 


11.74 


1.00 


1 


,243±0.018 


-0.053+0.012 






2MASS J05021 121-6759295 


11.97 


0.96 


1 


,197±0.026 


-0.045+0.024 






2MASS J05021 137-6758401 


11.98 


1.02 


1 


.199+0.023 


-0.064+0.016 






2MASS J06002748-6342222 


9.60 


1.80 


1 


.075+0.012 


0.437+0.005 


+,c 


NGC 2162 


2MASS J06003 156-6342581 


11.64 


1.03 


1 


.202+0.011 


-0.043+0.011 


+ 




2MASS J060033 1 6-6342 13 1 


12.24 


0.99 


1 


,197+0.016 


-0.127+0.017 


+ 




2MASS J06003869-6341393 


12.26 


0.99 


1 


.184+0.020 


-0.038+0.017 






2MASS J05563892-7258155 


11.77 


1.03 


1 


,227+0.012 


-0.091+0.013 




NGC 2173 


zMAs!) JUjj / 566 /- 1 Ijalyy 


IZ.Ul 


a nc 


1 

1 


1 Ct~ i A A 1 A 

. 166+U.U14 


A t\AH i A AT7 

-0.04/+U.U1 / 


+ 




2MASS J05570233-7257449 


12.13 


1.04 


1 


.159+0.016 


-0.008+0.019 






2MASS J05575784-7257548 


12.18 


1.03 


1 


.203+0.016 


-0.048+0.019 


+ 




2MASS J05563368-7257402 


12.43 


1.00 


1 


.146+0.016 


0.010+0.027 






2MASS J05581 142-7258328 


12.45 


1.04 


1 


.154+0.019 


-0.021+0.023 


+ 




2MASS J05583257-7258499 


12.48 


0.92 


1 


.136+0.015 


-0.070+0.022 






2MASS J05572334-7256006 


12.56 


1.01 


1 


.227+0.019 


-0.016+0.027 






2MASS J0556576 1-7254403 


12.84 


1.01 


1 


.169+0.020 


-0.099+0.031 







Notes. (1) 2MASS ca talogue star name, (2) extinction corrected A"-band magnitude and (3) (J-K) colour from the 2MASS Point Source Catalogue 
(Skrutskie et al. 2006), (4) K-band D C o index value, (5) //-band C2 index value, (6) Notes on individual stars: "C" - a carbon-rich stars, "+" - 
the star was used for the integ rated spectrum of the cluster, ( 7) globular cluste r, next to which the star was observed. The Deo index is defined in 
iMarmol-Oueralto et ail d2008f) and the C 2 index is taken from lMarastonl d2005l) . 



is consistent with the model predictions. We note that we were 
unable to measure the C2 index in the other two old and metal 
poor clusters in our sample, due to the very low quality of their 
//-band spectra. 

The group of the three intermediate age and more metal rich 
clusters is best compared with the SSP model with [Z/H] = 
-0.33 (blue, dotted line). The C2 index measurement for 
NGC 2173 (2Gyr, red solid symbol) is consistent with the 
model prediction of increasing C2 index towards younger ages. 
The youngest globular clusters at ~1.3Gyr (NGC 1806 and 
NGC 2162) follow this trend, although their C2 index is signif- 
icantly larger than the model prediction. The latter observation 
supports our findings based on the /iT-band Deo index in Paper I. 
There, we argued that the reason for the discrepancy between 
the CO index predictions of the models and the data (see Fig. |4] 
top panel) is due to the presence of carbon-rich stars which in- 
fluence the spectrophotometric properties of stellar populations 



at intermediate ages. The empirical calibration of the models is 
b ased on spectra of the Milky Way carbon stars from the library 
of lLancon & Mouhcind (120021) . combined with photometric cal- 
ibration based on the LMC globular clusters, which led to in- 
consistent results. The presence of LMC carbon-rich stars in the 
spectra leads to a decrease of the Deo index, but to an increase 
of the C2 index. To test this claim, we used the spectrum of NGC 
2173, which is the "oldest" intermediate-age globular cluster 
and therefore the one wi th the least contribution from carbon- 
rich stars in our sample dMucciarelli et al.l l2006). as a baseline 
for a simplistic stellar population modelling test. Increasing the 
contribution from TP-AGB stars in this cluster would mimic a 
younger age, such as seen in NGC 1806 and NGC 2162. In this 
way we are able to test the effects of different carbon- to oxygen- 
rich stars ratios in a given stellar population on its spectral fea- 
tures. When we added 40% to 70% of the averaged Milky Way 
stellar spectrum contained in the 3rd bin of lLancon & Mouhcind 
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Fi g. 4. Compariso n between C2 index values and the models 
of Marasto nl OOOl . The ages of the clusters are the mean ones 
from Table[T]with error bars reflecting the r.m.s. For guidance, in 
the top panel we reproduce the /if -band Dco index from Paper I 
(their Fig. 9). In both panels with black arrows we indicated 
how the spectrum of NGC2173 changes when adding differ- 
ent fractions of LMC carbon star light to it. The red arrows in- 
dicate the changes in the same spectrum, but adding different 
fractions of the averaged M ilky Way carbon star spectrum of 
lLancon & Mouhcind(l2002l) . See text for a detailed description. 



d2002l) to the spectrum of NGC 2173, the resulting C 2 and D C o 
indices increased, as predicted by the models for the younger 
clusters (red arrows in Fig. 4). When performing the same test, 
but adding the spectrum of a carbon star in the LMC (in this case 
the star 2MASS J06002748-6342222 in NGC 2162) to the clus- 
ter, the resultant Dco index decreased (black arrows in the same 
figure) and fitted the observations of the younger clusters better. 
Contrary to the Deo index dependence on C-star content, the re- 
sultant C2 index increased consistent with the model predictions 
regardless if the C-star is taken from the Milky Way or the LMC. 

The results from the above tests depend on the actual car- 
bon star spectrum that is used to perform the test. Carbon 
stars exhibit large variations in their properties for a given 
environment, as we discu ss further in Sect. |5] (see also e.g. 
lLancon & Mouhcmell2002h . This caveat illustrates the intrinsic 
problem of stochastic flu ctuations of AG B stars in globular clus- 
ters (see e.g. review bv lLanconll201 ll) . AGB stars are rare in 
stellar populations, with only one in a population of 10 4 stars 
in total. Nevertheless, a single star c an produce up to 80% of 
the stellar population's near-IR light (lMarastonll2005h . Here we 



see that it can also play a very important role in the integrated 
spectral properties of the hosting stellar population. This will 
inevitably lead to large or even significant deviations between 
the model predictions and observations of individual clusters. 
Various authors have made extensive studies of the influence 
of such stochast ic fluctuations for the integrated colours of stel- 
lar clusters (e.g. | Piskunov eta l. 2009; Popescu & Hansonll2010t 
iFouesneau & Lanconl2010l) . A similar study for the spectral fea- 
tures in integrated spectra of globular clusters would be benefi- 
cial for the interpretation and judgment of our results. 

5. Carbon-rich AGB stars in the LMC and the Milky 
Way 

A possible explanation for the differences between models and 
observations, described in the previous section, may be the dif- 
ferent metallicities or carbon and oxygen abundances of the sam- 
ple stars. In more metal p oor stars (i.e. th e LMC stars), the C/O 
ratio is higher on average. livlatsuura et al.l(l2002l) argue for a sys- 
tematically larger C/O ratio in LMC carbon stars compared to 
Galactic ones. Models of carbon-rich stars with high er C/O ratio 
have weaker CO features , but stronger C2 features dLoidl et al.l 
l200lHAringeretal.l2009h . We note that a variation in the molec- 
ular bands of carbon stars during their pulsation cycles is not ex- 
pected to be strong dLoidl et al.l 2001b . In Paper I we supported 
this scenario, i.e. the difference in metallicity to be the cause of 
the discrepancy, by showing that stars in the LMC and the Milky 
Way with similar (J - K) colours can hav e s ubstantially differen t 
D co index values, as lCohen et all dl98ll) and lFrogel etail (Il990h 
have shown earlier as well. 

Based on the limited stellar samples discussed there, we 
speculate that the LMC carbon-rich stars have weaker Dco in- 
dices when compared to carbon stars in the Milky Way. In 
Fig. [5] (left pa nel) we repeat the di agram and add stars from the 
SpeX library dRavner et al.l 120091 . This library contains spec- 
tra of mostly near-solar metallicity K and M stars with lumi- 
nosity classes between I and V, carbon-rich stars and S-stars 
(for which C/O^l). For half of the carbon-rich stars (four) we 
could find iron abun dance estimates in the PASTEL database 
(ISoubiran et al.l20loh . ranging from -0.3 to 0.2 dex. Their mean 
iron abundance is slightly sub-solar, which is on average more 
metal-rich than the sample of LMC carbon stars with [Fe/H] 
-0.4. The carbon-rich stars from the SpeX sample are located 
between the two relations, found in Paper I, for LMC (solid 
lin e) and Milky Way (dashed li ne) carbon stars from the library 
of lLancon & Mouhcinel d2002l) . thus smearing any clear rela- 
tions. Moreover, if instead of the averaged Milky Way c arbon 
star spectra, the individual stars of lLancon & Woodl (120001 were 
plotted, then the dispersion would have been even more promi- 
nent. However, we do not observe Dco in any of the Milky Way 
stars as weak as in three of the LMC stars. Surprisingly, on the 
(J — K)- C2 plot (Fig. [5] right panel) there is no obvious differ- 
ence between the stars in the two galaxies. Assuming the lower 
metallicity of the LM C and the higher C/O ratio, as argued by 
iMatsuura et al.l (120021) . the C2 index is expected to be stronger 
in the LMC than in the Milky Way carbon stars. Instead we find 
the opposite: that two of the Galactic stars have a stronger Co 
index than the LMC stars. This could be due to the C/O ratio 
changing at every dredge-up episode, thus leading to a disper- 
sion in the Co index larger than the metallicity effect. Based on 
the three stellar samples explored here, we cannot confirm that 
there is a marked difference between the Milky Way and LMC 
carbon stars in terms of the C/O ratio. To reach conclusive re- 
sults about the dependence of the C2 index on the C/O ratio, 
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different element abundances, and the metallicity of the stars, 
detailed abundance measurements for a larger sample of LMC 
and Milky Way giant stars are necessary. 

Fig.[6]shows the behaviour of the A'-band CO versus the H- 
band C2 features in stellar spectra. The complete sample of stars 
displays a large spread of Deo index values, while the C2 in- 
dex is significantly stronger than zero only in carbon-rich stars. 
There are also a number of M-type stars (M6-9III), whose C2 
index is larger than zero but is still systematically smaller than 
the value of Co in C-type stars. In these stars the re is significant 
H2O absorption in the ran ge 1.7 - 2.05 //m (e.g. lMatsuura et all 
119991; iRavner et al.l l2009t) that leads to a steep decrease of the 
continuum at the location where the C2 feature is located. Thus, 
when measuring the C2 index in these M giants we obtained a 
non-zero index value, while in reality there is no C2 absorption. 

Based on this figure, there is no clear separation between 
LMC and Milky Way carbon-stars, although the LMC stars have 
on average lower Dco indices, which is consistent with their 
lower metallicity. However, the same LMC stars do not show a 
stronger C2 index, which would be expected as discussed earlier 
in the text. 

Th ere is one C-type star from the catalogue of IRavner et al.l 
d2009l) that shows Deo an d C2 indices and a (/ - K) colour con- 
sistent with K and M-type stars (see Fig. [Hi. T his is HP 76846, 
classifi ed as C-R2 dKeenan & Newsom l2000h . As lTanaka et alJ 
J2007h have shown, the C2 feature is intrinsically weak in early 
hotter C-R stars. 

In Fig. |6l for comparison, we show the measured Deo an d 
C2 indices from the integrated LMC globular cluster spectra 
(coloured triangles). The two clusters that are dominated by car- 
bon stars, NGC 1806 and NGC2162, follow the stellar trends of 
higher C2 index values. Such a near-IR spectroscopic index di- 
agnostic plot can be very helpful in revealing the age of a given 
integrated stellar population, being a stellar cluster or a galaxy. 
As shown above, the C2 feature is present and its strength is 
significan t only in carbon star s and stellar populations harbour- 
ing such. lLancon et all (1 19991) have already proposed a similar 
C2 index, that can be used to detect the presence of intermedi- 
ate age stellar populations. Even when smeared out by the line- 
of-sight velocity distribution of stars in galaxies of the order of 
400 km s" 1 , it will remain identifiable and measurable. By ob- 
serving an increased strength of the C2 index, one can conclude 
that the object is in the carbon-star phase and thus date it to an 
intermediate age of 1-2 Gyr. 

6. Concluding remarks 

In this paper we concluded our pilot project of providing an em- 
pirical near-IR library of integrated spectra of globular clusters 
in the LM C by adding J- a nd //- band spectra to the /f-band 
ones from iLvubenova et al.l (1201 Oi Paper I). We provided full 
near-IR spectral coverage (as observable from the ground) for 
four globular clusters: the old and metal poor cluster NGC 20 19, 
and the intermediate age and approximately half-solar metallic- 
ity clusters NGC 1806, NGC 2162, and NGC 2173. For the old 
and metal poor clusters NGC 1754 and NGC 2005 the //-band 
part is missing due to bad observing conditions and unaccept- 
able sky emission subtraction residuals. Using our sample of 
globular clusters we tested predictions of current evolutionary 
population synthesis models in the near-IR and discussed some 
of the problems that arise from this comparison. Although still 
influenced by small number statistics, our spectra show that the 
near-IR spectral domain is a useful source of information when 
applied to the study of integrated stellar populations. 
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Fig. 6. Dco vs - C2 index for stellar and globular cluster spectra. 
Origin of spectra, symbols and colours as in Fig. [5] Error bars 
are about the size or smaller than the symbols, thus they are not 
explicitly shown. Coloured triangles show the measurements for 
the integrated spectra of LMC globular clusters. 

One of the main goals of our project was to provide a small 
empirical library of near-IR spectra of globular clusters to serve 
as a test bench for current and future stellar population models. 
In Fig. [2] we presented a comparison between the overall near- 
IR spectral energy dist ribution of sampl e clusters and the stellar 
population models of iMarastonl (I2005I) . The models agree re- 
markably well with the data in terms of the spectral shape and 
the most prominent absorption features. Especially, the predicted 
distinct signatures of the TP-AGB phase in intermediate age (1- 
2 Gyr) stellar populations are nicely reproduced. 

However, the agreement is not so good when exploring indi- 
vidual spectral features in more detail, where more subtle differ- 
ences in absorption strength can be quantified. Here, as well as in 
Paper I, we concentrated our studies to the strongest absorption 
features in the near-IR, namely C2 at 1 .77 fim and 12 CO (2-0) 
at 2.29 yum, due to the limited spectral resolution of the current 
models. We described their behaviour using two indices: C2 and 
Dco, respectively. The first one reflects the presence of carbon- 
rich AGB stars in the stellar population, while the second is typ- 
ical for both carbon- and oxygen- rich AGB and RGB stars. We 
refer the reader to Sect. 2] where we discussed in detail the ob- 
se rved differences b etween our observations and the EPS models 
of lMarastonl d2005l) . In Sect. [5] we explored some possible rea- 
sons for these, based on stellar libraries with Milky Way and 
LMC carbon-rich stars (see Fig. |5). These stars can contribute 
up to 80% to the total /if -band light of a globular cluster stellar 
population, thus their proper inclusion in models is crucial. Our 
results were inconclusive about the proposed difference in the 
C/O ratio between stars in the Galaxy and the LMC. To confirm 
or reject this, accurate abundance measurements of iron, oxygen 
and carbon elements are necessary. Therefore, there is a clear 
need for larger spectral libraries of carbon-rich stars with a range 
of metallicities before we are able to reproduce their proper con- 
tribution to the stellar population models. 

The occurrence of a strong C2 index in the near-IR spectra 
of globular clusters and hence in galaxies is an indication that 
they are in the carbon star phase, which dates the hosting stellar 
popula tion to 1 - 2 Gyr, as previously suggested by lLancon et al.l 
dl999l) . In contrast, the Dco index is much less straightforward 
to interpret, since the magnitude of this index is driven by a 
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Fig. 5. Deo (left) an d Co (right panel) indices versus (J - K) colour for stars in the Milky Way and the LMC. The oxygen- and 
carbon-rich LMC stars (filled black and red dots, respectively) are from our own SINFONI observations, the K (green o pen squares), 
M (bla ck diamonds), C (red asterisks), and S (red asterisks with square around) stars are taken from the SpeX library of iRavner et al.l 
d2009l) . The C-rich spectra in the Milky Way (orange crosses) are taken from the library of lLancon & Mouhcinel (12002 ). Typical 
error bars are s hown in the upper right corner of the left panel. In the left panel, the solid line denotes the relation for LMC carbon- 
rich stars fromlLvubenova et al.l (1201 Ol Eq. 3), the dashed line (Eq. 5 in the same paper) represents the relation for the Milky Way 
spectra of lLancon & Mouhcinel (120021) . 



complex combination of metallicity, surface gravity (luminos- 
ity), and effective temperature, as illustrated by Fig. In Fig. [6] 
we propose a diagnostic plot, based only on these two near-IR 
indices. Both of them are sufficiently broad and strong, so they 
can be easily measured in galaxy spectra, even when smeared out 
by the line-of-sight velocity distribution of the stars. For galaxies 
up to a redshift of 0.007, e.g. Virgo and Fornax clusters, the C2 
feature is in the rest-frame //-band. For redshifts between 0.007 
and 0.13 the feature is hidden by the atmospheric cutoff, but it 
remains accessible from space. 

During the next decades new facilities, such as the JWST 
and the E-ELT, will widen the discovery space in extragalactic 
research. Their enhanced sensitivity will be in the infrared por- 
tion of the spectrum. However, the stellar population analysis 
methods for the integrated light are not yet fully developed in the 
near-IR wavelength regime. The methods discussed here offer a 
step towards a better understanding of galaxy evolution and ob- 
servations are already feasible with current facilities. For exam- 
ple, current adaptive optics systems, working mostly in the near- 
IR, offer a possibility to resolve the innermost parts of nearby 
galaxies. New information about the stellar population proper- 
ties at these spatial scales can give us important clues about how 
galaxies have formed and evolved. 
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